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Abstract 



We report high precision transit photometry of GJ1214b in JHK S bands taken simultaneously with the 
SIRIUS camera on the IRSF 1.4 m telescope at Sutherland, South Africa. Our MCMC analyses show 
that the observed planet-to-star radius ratios in JHK S bands are R p /R s j = 0.11833 ± 0.00077, i? p /i? s .H = 
0.11522 ±0.00079, R p /R B , Ks = 0.11459 ±0.00099, respectively. The radius ratios are well consistent with 
■ the previous studies by Bean et al. (2011) within ler, while our result in K s band is shallower than and 

Qh inconsistent at 4<r level with the previous measurements in the same band by Croll et al. (2011). We have no 

q . good explanation for this discrepancy at this point. Our overall results support a flat transmission spectrum 

in the observed bands, which can be explained by a water-dominated atmosphere or an atmosphere with 
r/3 , extensive high-altitude clouds or haze. To solve the discrepancy of the radius ratios and to discriminate a 

definitive atmosphere model for GJ1214b in the future, further transit observations around K s band would 
be especially important. 
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1. Introduction to as transmission spectroscopy, is useful to learn exoplan- 

etary atmospheres (e.g., Charbonncau et al. 2002; Winn 

Among several hundreds of discovered extrasolar plan- et al. 2004; Narita et al. 2005; Redfield et al. 2008; Snellen 

ets, transiting exoplanets can give us the most informa- et al. 2008; Swain et al. 2008; Desert et al. 2009; Sing et al. 

' tion about properties of planets. Transiting exoplan- 2011). 

ets are unique in that both the true mass and the ra- Recent discoveries of transiting planets with smaller 

dius can be determined by measurements of radial veloc- radii (i.e., a few Earth radii) around M dwarfs have ex- 

ities and transit depths. The mass-radius relation en- panded targets of this methodology down to terrestrial 

■ ables us to infer the internal structure and bulk com- exoplanets. Transiting planets around M dwarfs are espe- 

j_h positions of transiting exoplanets (e.g., Valencia et al. cially favorable for this kind of studies since transit depths 

5^ ■ 2007). One can also measure the angle between the stel- become relatively deeper due to smaller host stars' radii, 

lar spin axis and the planetary orbital axis in the sky In this sense, GJ1214b discovered by Charbonneau et al. 

projection via the Rossiter-McLaughlin effect, which pro- (2009) is currently the most interesting target. Numbers 

vides insights of planetary migration histories (e.g., Winn of previous observers reported transit depths of GJ1214b 

et al. 2005; Narita et al. 2007; Hebrard et al. 2008; Narita in various wavelength bands (e.g., Bean et al. 2010; Bean 

et al. 2009; Winn et al. 2009; Triaud 2011; Hirano et al. et al. 2011; Croll et al. 2011; Desert et al. 2011; Carter 

2011; Albrecht et al. 2012). Another great merit of et al. 2011; Berta et al. 2011; Bcrta et al. 2012; de Mooij 

transiting planets is that they also provide opportunities et al. 2012). The previous observations have revealed 

to probe atmospheric compositions by transit photome- that there are two possible atmospheric models to explain 

try or spectroscopy. Theoretical predictions have shown the observed transit depths: a water-dominated (higher 

that transit depths depend on wavelengths (e.g., Seager mean molecular weight / lower scale height) atmosphere 

& Sasselov 2000; Millcr-Ricci & Fortncy 2010; Howe & and a hydrogen-dominated (lower mean molecular weight 

Burrows 2012), reflecting atmospheric compositions and / higher scale height) atmosphere. Among the previous 

environments (e.g., haze/clouds). Previous space-based studies, supporting evidences for the hydrogen-dominated 

and ground-based observations for transiting hot Jupiters atmosphere are deeper transit depths observed in K s band 

have indeed revealed that this methodology, often referred (Croll et al. 2011; de Mooij et al. 2012) and in g band (de 
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Mooij et al. 2012). The rest of observations have shown a 
fairly flat transmission spectrum for GJ1214b through vis- 
ible and infrared wavelengths, which indicates the water- 
dominated atmosphere or the hydrogen-dominated but 
with high- altitude clouds/haze atmosphere. Thus the na- 
ture of this planetary atmosphere is still an open question. 

One known problem of GJ1214 as a target of transit 
photometry is possible stellar variability and starspots, 
which may cause small time-dependent variations in tran- 
sit depths. Some of the previous studies estimated that 
such small variations would correspond to the ratio of 
radii of the planet and star R p /R s ~ 0.001 (Carter et al. 
2011; Bcrta et al. 2011; dc Mooij et al. 2012), which are 
indeed comparable to or larger than the nominal observa- 
tional uncertainties reported in the previous publications. 
This fact makes it difficult to compare the all observed 
transit depths in different epochs. For this reason, si- 
multaneous multi-band transit photometry is the most 
useful way to discriminate possible atmosphere models. 
We thus focus on JHK S simultaneous transit photometry 
since the largest difference in transit depths is theoreti- 
cally predicted between JH bands and K s band. In this 
sense, our motivation is the same as that by Croll et al. 
(2011), who tried nearly simultaneous J and K s band tran- 
sit photometry by rapidly switching the two bandpass fil- 
ters. Interestingly, Croll et al. (2011) claimed that K s 
band transit depth is significantly (5ct) deeper than that 
in J band. With sufficient precision, we can independently 
confirm or refute the claim by Croll et al. (2011) via JHK S 
simultaneous transit photometry. 

In this paper, we present a result of such an observa- 
tion using the SIRIUS camera onboard the IRSF 1.4 m 
telescope. This is indeed the first high precision tran- 
sit observation using IRSF/SIRIUS. We summarize target 
properties and observation methods in section 2, and de- 
scribe our analysis methods in section 3. We report our 
analysis results in section 4, and discuss its meaning in 
section 5. Finally, we summarize this paper in section 6. 

2. IRSF SIRIUS Observation for GJ1214 

2.1. Target Properties 

The host star GJ1214 is an M4.5 type star at 13 pc 
away from the Sun, with a mass of 0.157±0.019 Mq 
and a radius of 0.2110 ± 0.0097 R Q ( The star is known 
to show stellar variability at ~2% level with a period of 
~52.3 d (see e.g., Berta et al. 2011) and also known to 
have starspots with surface coverage variability at several 
percent level (Carter et al. 2011). 

The orbiting planet GJ1214b is categorized as a super- 
Earth with a mass of 6.55 ± 0.98 Mq and a radius of 
2.68 ±0.13 R B (Charbonncau et al. 2009). The planet 
orbits around GJ1214 at a period of ~1.58 d and a semi- 
major axis of ^0.0146 AU, where is slightly closer to the 
host star than the star's habitable zone (ohz ~ 0.06 AU). 
Previous investigation of transit timing variations (TTV) 
has shown no evidence of large TTV in this system (Carter 
et al. 2011). 



Table 1. Photometric transit light curves of GJ1214b taken 
by IRSF/SIRIUS 

"Tunc [BJDtdb] Value Error 
J band 



2455788.28330 
2455788.28396 
2455788.28462 
2455788.28527 
2455788.28593 



0.99979 0.00124 

1.00011 0.00124 

1.00125 0.00124 

1.00137 0.00125 

1.00015 0.00125 



All data arc presented in the electric table. 



2.2. Observation Setup 

We observed a full transit of GJ1214b during 18:43- 
21:00 of UT 2011 August 14 with Simultaneous Infrared 
Imager for Unbiased Survey (SIRIUS: Nagayama et al. 
2003) on the Infrared Survey Facility (IRSF) 1.4 m tele- 
scope. The IRSF and the SIRIUS camera were con- 
structed and has been operated by Nagoya University, 
SAAO (South African Astronomical Observatory) at 
Sutherland, South Africa, and National Astronomical 
Observatory of Japan. The SIRIUS camera utilizes two 
dichroic filters and three 1024x1024 HgCdTe detectors, 
which can obtain J (1.25/im ± 0.085/im), H (1.63/xm± 
0.15/im), K s (2.14/im ± 0.16/^m) band images simultane- 
ously, with a square field of view of 7.7' on a side and 
a pixel scale of 0.45". To achieve higher photometric 
precision, we defocused the telescope so that stars have 
doughnut-like point spread function (PSF). The typical 
size of the PSF was ~ 13 pixels in radius. We carefully 
set GJ1214 and bright comparison stars away from bad 
pixels on the detectors. The exposure time was 49 s and 
the typical dead time was 7.6 s (duty cycle of 86.6%). As 
a result, we obtained 145 frames during the above time 
and use them for subsequent analyses. 

2.3. Position- Locking Software 

For high precision infrared photometry, it is known 
that target positions on the detectors should be stared 
rather than dithered (see e.g., Croll et al. 2011). Since 
the IRSF does not equip an auto-guider, target positions 
tend to move over tens of pixels during observations and 
it was the biggest problem for precise transit observations 
with the IRSF. To avoid such changes of target positions, 
we have developed and installed a customized position- 
locking software which gives appropriate feedback to the 
telescope when the ccntroid position of the target star 
moves over 1 pixel on the detectors. By using the position- 
locking software, GJ1214's centroid positions on the three 
detectors are kept within an rms of about 2 pixels in both 
x and y directions. 

3. Data Reduction and Light Curve Analysis 

Primary data reduction is carried out with a pipeline 
for the SIRIUS 1 , including a correction for non-linearity, 



http:/ /www. z. phys.nagoya-u.ac.jp/~nakajima/sirius/software/software. html 
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dark subtraction, and flat fielding. We note that the de- 
tectors of SIRIUS saturate at 32000 ADU (J band detec- 
tor), 33000 ADU (H band detector), and 25000 ADU (K s 
band detector) respectively under the Correlated Double- 
Sampling (CDS) readout mode, and keep a good linearity 
(<1%) up to -10000 ADU. The above non-linearity cor- 
rection enables us to work up to —25000 ADU with <1% 
linearity. The correction is necessary to reduce possible 
systematic errors stemmed from non-linearity in fractional 
light curves down to <0.1%. This is especially important 
for H and K s bands where sky background levels arc high. 

Aperture photometry is followed using a customized 
aperture photometry pipeline (Fukui ct al. 2011), with 
a constant aperture-size mode where a same aperture size 
is applied for all images. Sky background levels on the ob- 
serving night were typically ~ 1000 ADU, - 7000 ADU, 
and —6000 ADU, for JHK S bands, respectively. Maximum 
counts in the frames (including sky background) are below 
25000 ADU, where the non-linearity correction can work. 

Time stamps of observations are recorded in the FITS 
headers in units of Modified Julian Day (MJD) based 
on Coordinated Universal Time (UTC). We convert the 
time system to Barycentric Julian Day (BJD) based on 
Barycentric Dynamical Time (TDB) using the code by 
Eastman et al. (2010). 

We then select appropriate comparison stars and aper- 
ture radii for each band so that light curves have the small- 
est rms at out-of-transit (OOT) phase as follows. There 
are three possible comparison stars within the field of view 
of SIRIUS: one is almost the same brightness as GJ1214 
and the others are 1-2 magnitude fainter than GJ1214. 
First, we create fractional light curves by dividing fluxes 
of GJ1214 by sum of all combinations of comparison stars' 
fluxes with changing the aperture radius from 13 pixels to 
30 pixels in increments of 1 pixel. All the fractional light 
curves fobs exhibit systematic trends (0.1% level and ap- 
parently linear) at OOT phase. The systematic trends 
could arise due to slow variability in the brightness of 
GJ1214 itself or comparison stars, changing airmass, or 
position changes of the stars on the detectors, and so on. 
We then correct the systematic trends by the following 
expression, F cor = F ODS x 10 _0 - 4Amcor , by assuming the 
correction factor Am cor = ko + k t t + k z z + k x dx + k y dy, 
where t is the time, z is the airmass, dx and dy are the 
differences of centroid positions in x and y directions, and 
ko,k tl k z ,k x ,k y are free coefficients. We eliminate outliers 
(2, 10, 1 frames for JHK S bands, respectively) with large 
centroid position changes over 10 pixels, and determine 
optimal coefficients using OOT data by the AMOEBA al- 
gorithm (Press et al. 1992) so that the corrected OOT data 
have the smallest rms. In this process, we also rescale the 
photometric errors of the data so that reduced \ 2 become 
unity. We consequently choose to use only the brightest 
comparison star (the same one for all bands) and aper- 
ture radii of 15, 16, 15 pixels for JHK S bands, which give 
the smallest rms of the data as 0.00124, 0.00125, 0.00155 
for JHK S bands, respectively. We note that including 
the other two comparison stars or using different aper- 
ture radii have less impact on our subsequent conclusion, 
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but with slightly larger errors. The data of the corrected 
transit light curves are presented in table 1. 

For modeling the corrected transit light curves, we fol- 
low the procedures that were adopted in previous stud- 
ies in litelature (Bean et al. 2011; Croll et al. 2011; dc 
Mooij et al. 2012; Berta et al. 2012) so as to compare 
our transit depths with previous ones. We fix the or- 
bital inclination i to 88.94° and the orbital distance in 
units of the stellar radius a/R s to 14.9749, which were 
determined by Bean et al. (2010). We also fix the or- 
bital period of GJ1214b (P = 1.58040481 d) and the mid- 
transit time (T C)0 = 2454966.525123 BJD T db), determined 
by Bean et al. (2011). This assumption is justified by 
the fact that there is no evidence of large TTV (Carter 
et al. 2011). We use the quadratic limb-darkening law, 
I(fi) = 1 — ui(l — yu) — 1*2(1 — /i) 2 , where I is the inten- 
sity and [i is the cosine of the angle between the line 
of sight and the line from the stellar center to the po- 
sition of the stellar surface. We adopt empirical quadratic 
limb-darkening coefficients (mi,j = 0.088, it2,j = 0.404, 
ui )H = 0.076, u 2 ,H = 0.407, ui, Ks '= 0.048, u 2 ,k s = 0.350) for 
JHK S bands given by Claret & Bloemen (2011), assum- 
ing stellar effective temperature T e g = 3000 K and log of 
the stellar surface gravity log g = 5.0. We note that these 
assumptions on T e g and logg are the same as Croll et al. 
(2011) for J and K s bands, although Croll ct al. (2011) 
used the non-linear limb-darkening law. 

We fit the JHK S transit light curves simultaneously us- 
ing the analytic formula given by Ohta et al. (2009), which 
is equivalent with Mandel & Agol (2002) when using the 
quadratic limb-darkening law. The four free parameters 
are the difference of the observed mid-transit time from 
the predicted mid-transit time AT C for this epoch (epoch 
E = 520 from T c .o) and the radius ratios of the planet 
and star R p /R s for JHK S bands. In addition to the main 
analysis, we also test some other cases with different con- 
ditions as follows: (1) T c for each band to be free, (2) 1*2 
for each band to be free, and (3) i and a/R s to be free, 
to assess the validness of assumptions of fixed parame- 
ters above. We employ the Markov Chain Monte Carlo 
(MCMC) method to estimate values and uncertainties of 
the free parameters, following the analysis by Croll et al. 
(2011). We create 5 chains of 5,000,000 points, then trim 
the first 100,000 points from each chain. The acceptance 
ratios of jumping for the chains are set to about 25%, and 
we check that they pass the Gelman & Rubin (1992) test. 
We define la statistical errors by the range of parameters 
between 15.87% and 84.13% of merged posterior distribu- 
tions. 

4. Results 

The mean values and la statistical errors based on 
the MCMC analyses are summarized in table 2. We 
note that the errors do not include systematic errors dis- 
cussed in the next section. In the main analysis, we 
find AT C = 0.000104 ± 0.000067 s, which is consistent 
with the predicted mid-transit time within 2a, validat- 
ing the assumption to use the fixed period determined 



JHK S Simultaneous Transit Photometry of GJ1214b 



4 N. Narita et al. [Vol. , 



Table 2. Mean values and errors of the parameters based on the MCMC analyses. 



Parameter 


Value 


Error^ 


Value 


Error 


Value 


Error 


Value 


Error 




main 


case 


(1) 


case 


(2) 


case 


(3) 


-R p /i? Sj j 


1 1833 


+0 00077 


0.11833 


±0.00077 


0.11867 


±0.00095 


0.11800 


±0.00083 


Rpj i? s ,H 


0.11522 


±0.00079 


0.11519 


±0.00079 


0.11668 


±0.00096 


0.11492 


±0.00085 


Rpl fis,Ks 


0.11459 


±0.00099 


0.11463 


±0.00099 


0.11428 


±0.00121 


0.11427 


±0.00101 


AT 0>aU [s] 


0.000104 


±0.000067 


- 


- 


0.000107 


±0.000066 


0.000102 


±0.000069 


AT CiJ [s] 






0.000111 


±0.000111 


- 


- 


- 


- 


AT CiH [s] 






0.000005 


1 i~\ AAA1 "1 A 

±0.000114 










AT c ,ks [s] 






0.000244 


±0.000134 










U2, J 


0.404 


fixed 


0.404 


fixed 


0.345 


±0.095 


0.404 


fixed 




0.407 


fixed 


0.407 


fixed 


0.145 


±0.103 


0.407 


fixed 


«2,Ks 


0.350 


fixed 


0.350 


fixed 


0.403 


±0.125 


0.350 


fixed 


Ml 


88.94 


fixed 


88.94 


fixed 


88.94 


fixed 


89.13 


+0.53 
-0.44 


a/R s 


14.9749 


fixed 


14.9749 


fixed 


14.9749 


fixed 


15.0468 


+0.3413 
-0.4847 



t The presented errors do not include systematic errors due to the fixed parameters (sec text). 
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Fig. 1. Transit light curves with the best-fit models in J (top 
panel), H (middle panel), K s (bottom panel) bands, observed 
simultaneously by the IRSF/SIRIUS. Each lower panel shows 
residuals between the observed values and the best-fit models. 
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Fig. 2. MCMC histograms of the free parameters. The bin- 
ning size is set to a hundredth part of the plotting region. The 
dashed vertical line in each panel indicates the mean value. 
The thick zones represent the la statistical errors shown in 
table 2. 

by Bean et al. (2011). The derived radius ratios are 
-R P /i? s ,j = 0.11833±0.00077, i? p /i? s ,H = 0.11522±0.00079, 
R P /R s ,Ks — 0.11459 ±0.00099, respectively. Comparisons 
of our result with the previous studies are discussed in the 
next section. Figure 1 plots the transit light curves with 
the models using the mean MCMC values, and figure 2 
shows histograms of the posterior distribution of the free 
parameters. 

One may wonder why our errors are comparable with 
Croll et al. (2011), even though we observed only one tran- 
sit with the IRSF 1.4m telescope and Croll et al. (2011) 
observed 3 transits with the CFHT 3.8m telescope. The 
reasons are that Croll et al. (2011) switched between J and 
K s band filters during the transits and the 3.8m telescope 
was too large for the brightness of GJ1214. Namely, Croll 
et al. (2011) employed 4 s exposure and the total duty 
cycle was only 22% for J and K s bands, which is signif- 
icantly lower than our duty cycle of 86.6%. Considering 
the above fact, it is reasonable that our errors are similar 
to those by Croll et al. (2011). 
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From the case (1), we find that AT C for JHK S bands are 
all consistent with zero within 2a and one another. Also, 
the result shows that fitting AT C for JHK S bands at once 
or separately has less impact on radius ratios. 

The case (2) illustrates a potential of systematic errors 
and underestimated errors stemmed from the fixed limb- 
darkening coefficients. The derived U2 for J and K s bands 
are well consistent with the empirical values by Claret & 
Bloemen (2011), and then the radius ratios in the same 
bands are also well accorded with the main result. While 
the derived «2 in H band is — 2.5er apart from the empir- 
ical value, and the radius ratio is also different from the 
main result by ~1.5<r. Although some systematic errors 
may be present, we consider that the assumption of the 
adopted limb-darkening coefficients is reasonable for the 
mean values, since the empirical values are in excellent 
agreement with the derived values in J and K s bands. For 
this reason, we believe that the radius ratio in H band 
using the empirical limb-darkening coefficients (the main 
case) is more reliable than that with free (the case 2). 
Another fact seen from the case (2) is that the errors of 
the radius ratios are about 20% larger than the main case 
where the limb-darkening coefficients are fixed. This fact 
is often pointed out in previous studies (e.g., Southworth 
2008). Thus we note that the quoted errors in the radius 
ratios with the fixed limb-darkening coefficients may have 
20% additional errors. 

Finally, we confirm that the orbital inclination ({) and 
the orbital distance in units of the stellar radius (a/R s ) 
derived from our IRSF data are well consistent with those 
by Bean et al. (2010) via the analysis of case (3). We also 
confirm the derived radius ratios are consistent with the 
main result, however we notice that the errors are about 
10% larger than the main result. This fact also suggests 
that errors may be about 10% underestimated when using 
the fixed i and a/R s . 

5. Discussions 

5.1. Impacts of Systematic Errors 

In this section, we compare our radius ratios in JHK S 
bands with the previous studies and atmospheric models. 
Before doing that, we should care about possible system- 
atic errors. 

First, as assessed by the test cases (2) and (3), we note 
that our statistical errors may have additional 30% un- 
certainties due to the systematic errors stemmed from the 
fixed parameters (20% from fixing the limb-darkening co- 
efficients, and 10% from fixing the orbital inclination and 
the orbital distance). 

Second, it is also necessary to consider impacts of stellar 
flux variability and stellar spots on our results so as to 
compare radius ratios with previous studies measured in 
different epochs and wavelengths. Desert et al. (2011) 
presented an equation for estimating a difference of radius 
ratio A(i?.p/i? s ) caused by stellar flux variability due to 
unocculted spots (see equation 7 of Desert et al. 2011). 
The equation can be rewritten as 



A(R P /R S ) ~ 0.5 A/(A) (Rp/Rs), (1) 

where A/(A) is stellar flux variability at wavelength A. 
Berta et al. (2011) reported stellar flux variability of ~ 2% 
level in the MEarth observing bandpass (~ 780nm), and 
if following the assumption used in Croll et al. (2011) that 
GJ1214 (T off - 3000K) has spots 500 K cooler, then the 
stellar flux variabilities in JHK S bands can translate into 
- 1.5%,- 1.3%,- 1.0%, respectively (Croll et al. 2011). 
Using the equation above, possible variability of radius 
ratios in JHK S bands are - 0.75%, ~ 0.65%, - 0.5% of 
the observed radius ratios. Thus we should take into 
account possible systematic differences of A(i? p /i? Sj j) ~ 
0.00089, A(i? p /i? s , H ) - 0.00075, A{R p /R s<Ks ) ~ 0.00057, 
when comparing our radius ratios with other epochs' re- 
sults. 

Third, although we do not find any clear sign of spot 
occultation in our data, possible existence of unocculted 
spots is still a problem when comparing radius ratios 
with theoretical atmospheric models, since their impact 
depends not only on wavelength but also on spot prop- 
erties (coverage, temperature, and so on) and GJ1214's 
spot properties are still unclear, de Mooij et al. (2012) 
presented corrections for radius ratios assuming several 
cases of different spot coverage. According to their study, 
the corrections vary greatly between optical and infrared 
wavelength (A(R p /R s ) ~ 0.003 if spot coverage is 5%), but 
corrections in JHK S bands are very similar (see figure 7 of 
de Mooij et al. 2012). Thus we neglect the corrections for 
the current study since our data are only in JHK S bands, 
although we should take care when comparing our radius 
ratios with those in optical bands. 

5.2. Comparisons with the Previous Studies 

Figure 3 plots our observed radius ratios with those 
in the previous studies (Bean et al. 2011; Croll et al. 
2011; Desert et al. 2011; Berta et al. 2012; de Mooij et al. 
2012) and the two atmospheric models (i.e., hydrogen- 
dominated and water-dominated atmospheres) by Millcr- 
Ricci & Fortney (2010). Note that the values and errors of 
the previous studies are reported ones without any correc- 
tions, and the errors of our data are based on the MCMC 
analysis (30% additional errors are not included, since pre- 
vious observers did not take into account such systematic 
errors). 

In summary, our results support a featureless transmis- 
sion spectrum with a shallower (without a deeper) transit 
in K s band, which favors a water-dominated atmosphere 
rather than a cloudless hydrogen-dominated atmosphere 
by Miller- Ricci & Fortney (2010). A similar flat transmis- 
sion spectrum can be expected if the atmosphere is cov- 
ered by optically thick high-altitude clouds or haze (see 
e.g., Bean et al. 2011; Berta et al. 2012). 

A similar conclusion was also reported by Bean et al. 
(2011), who observed two transits of GJ1214b using the 
MMIRS instrument of the Magellan telescope by the 
multi-object spectroscopy approach. They obtained ra- 
dius ratios of GJ1214b in JHK bands (almost the same 
bandpass with our bands) as R p /R s j = 0.1158 ± 0.0024, 
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Fig. 3. Comparisons of our results with the previous ob- 
servational studies and the cloudless atmospheric models by 
Miller- Ricci & Fortncy (2010). The gray solid line represents 
the water-dominated (vapor) atmosphere model, and the gray 
dotted line does the hydrogen-dominated (Solar composition) 
model. The upper figure shows 1-2.5 fim wavelength region 
and the lower one does 0.6-4.5 /jm region. 

R P /R s ,n = 0.1 146 ± 0.0014, R p /R a>K = 0.1 158 ±0.0006, re- 
spectively. Thus our radius ratios in JHK S bands are all 
consistent with Bean et al. (2011) within a square-root of 
sum of squares of both la errors. Our radius ratios in 
J and H bands are also in agreement with Berta et al. 
(2012), who conducted spectro-photometry in the band- 
pass between 1.1 and 1.7 /im using the Wide Field Camera 
3 (WFC3) onboard the Hubble Space Telescope. 

On the other hand, our radius ratio in K s band appears 
inconsistent with the previous measurements in the same 
band by Croll et al. (2011) (R p /R s , Ks = 0.1199 ± 0.0008) 
and de Mooij et al. (2012) (i?p/i? s , Ks = 0.1189 ± 0.0015). 
Especially, our result (R p /R s . Ks = 6.1146 ± 0.0010) differs 
about 4cr (la here is a square-root of sum of squares of 
both la errors) from that by Croll et al. (2011). The 
discrepancy cannot be explained by stellar flux variability 
in K s band (A(i? p /i? s , Ks ) ~ 0.00057) nor 30% additional 
systematic errors due to the fixed parameters. Thus we 
do not confirm a deeper transit in K s band and we have 
no clear explanation for this discrepancy at this point in 
time. Further transit observations in K s band would be 
necessary to solve this enigma. 



5.3. Possible Origins of Water- Dominated Atmospheres 

The investigation of planetary atmospheres provides 
a crucial clue to learn the bulk composition of super- 
Earths like GJ1214b, which also gives important con- 
straints on formation processes of such planets. Only with 
its measured mass and radius, we are unable to distinguish 
whether GJ1214b is a rocky planet with a thick H/He at- 
mosphere or a water-dominated atmosphere (Rogers & 
Seager 2010; Nettelmann et al. 2011). Transmission spec- 
troscopy is thus an important key to discriminate the at- 
mospheric nature of transiting planets. As discussed ear- 
lier, our finding as well as some of the previous studies 
(e.g., Bean et al. 2011; Desert et al. 2011; Berta et al. 
2012) support a possibility of a water-dominated atmo- 
sphere for GJ1214b. Then, how can such atmospheres 
form? 

According to detailed modeling of the internal structure 
by Nettelmann et al. (2011), however, two-layer models 
(a water envelope on top of a rock core) which are con- 
sistent with GJ1214's age would have unreasonably high 
water-to-rock ratio. To fix this problem, a small amount of 
H/He, which accounts for several percent of the planet's 
total mass, must be incorporated in the water envelope 
(see Figure 4 of Nettelmann et al. 2011). In other words, 
GJ1214b is a Uranus/Neptune-like planet whose atmo- 
sphere is heavily polluted by water. 

Such properties of GJ1214b are qualitatively consistent 
with recent theories of planet formation. GJ1214b is or- 
biting close to its host star, in contrast to Uranus and 
Neptune. Its current location was too hot for icy build- 
ing blocks to exist in the proto-planetary disk where the 
planet was born. Therefore, provided GJ1214b contains a 
significant amount of water, it is reasonable to think that 
the planet has experienced orbital migration. 

One possible mechanism for migration is angular- 
momentum exchange with the proto-planetary disk (i.e., 
the type-I migration; Ward 1986). Recent iV-body simu- 
lations of planetary accretion with the effect of the type- 
I migration have demonstrated that water-rich proto- 
plancts of super-Earth-size can form in the vicinity of the 
central star (Ogihara & Ida 2009). As the proto-planets 
are embedded in the disk, they also capture the ambi- 
ent H/He disk gas naturally. A proto- planet with mass 
similar to that of GJ1214b can gain H/He atmosphere 
comparable in mass with that predicted by Nettelmann 
et al. (2011), although it depends on values of several pa- 
rameters (Ikoma & Hori 2012). Then, bombardments of 
smaller planetary embryos and/or other proto-planets rich 
in water can heavily pollute the H/Hc atmosphere by wa- 
ter during and after the migration. 

Another possibility is that a Uranus/Neptune-like 
planet that had formed beyond the snow line has migrated 
by gravitational scattering or by long-term interaction 
with outer gravitationally perturbing objects. Beyond the 
snow line, icy planctcsimals are sufficiently available for 
polluting H/He atmospheres. Indeed, the metallicity of 
the atmospheres of Uranus and Neptune are known to be 
significantly super-solar (Hubbard et al. 1995), which sug- 
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gest such pollutions. 

While both scenarios are to be quantitatively verified, 
the both scenarios suggest that there are other planets yet 
undetected in the GJ1214b system, although no TTV sign 
has been found (e.g., Carter et al. 2011). Future search for 
additional planets with RV measurements (possibly using 
brighter NIR region) is expected to verify the above sce- 
narios. Moreover, a measurement of the spin-orbit align- 
ment angle of GJ1214b via the Rossiter-McLaughlin effect 
would be very helpful in knowing which of the above two 
scenarios is true. 

6. Summary 

We conducted high precision JHK S simultaneous tran- 
sit photometry for GJ1214b using the SIRIUS camera on 
the IRSF at South Africa. It was the first high preci- 
sion transit observation using IRSF/SIRIUS. In this pro- 
cess, we have demonstrated that high precision NIR tran- 
sit photometry is possible even without an auto-guider by 
the position-locking software installed on the telescope. 
We have found a featureless transmission spectrum in the 
observed bands, which suggests a water-dominated at- 
mosphere or an atmosphere with extensive high-altitude 
clouds/haze. Although the observed radius ratios are well 
consistent with the results by Bean et al. (2011) and Berta 
et al. (2012), our result in K s band is unaccountably in- 
consistent with the previous studies in the same band 
by Croll et al. (2011). We have no good explanation 
for this discrepancy at this point. To solve the discrep- 
ancy and to distinguish a definitive atmosphere model for 
this planet, further transit observations around K s band 
would be especially important. In addition, it will be also 
interesting to observe transits in blue optical bands in 
order to ascertain the Rayleigh scattering in the atmo- 
sphere of GJ1214b. With such additional observations, 
we will be able to learn the true nature of the atmosphere 
of GJ1214b, and discriminating the true atmospheric na- 
ture would also give us useful insights on the formation 
process of this planet. 
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